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The role of the gluonic gg <-> ggg interactions in early 
thermalization in ultrarelativistic heavy-ion collisions 
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Abstract. We "quantify" the role of elastic as well as inelastic gg <-> ggg pQCD processes in kinetic 
equilibration within a pQCD inspired parton cascade. The contributions of different processes to kinetic 
equilibration are manifested by the transport collision rates. We find that in a central Au+Au collision 
at RHIC energy pQCD Bremstrahlung processes are much more efficient for momentum isotropization 
compared to elastic scatterings. For the parameters chosen the ratio of their transport collision rates 
amounts to 5 : 1. 

PACS. 05.60.-k Transport processes - 25.75.-q Relativistic heavy-ion collisions - 24.10.Lx Monte Carlo 
simulations 



1 Introduction 

It was speculated that a strongly coupled quark-gluon 
plasma (sQGP) pQ is formed in Au+Au collisions at RHIC. 
This strong coupling, or strong interaction, makes the 
QGP to a fluid with very small viscosity. However, how 
strong the coupling must be in order to generate a quasi- 
ideal fluid is an open question. 

Recently we have developed an on-shell parton cascade 
including elastic as well as inelastic gg <-> ggg pQCD pro- 
cesses to study the issue of thermalization 2 . Although 
the total cross section of the pQCD scatterings is a few mb, 
it is enough to drive the system into thermal equilibrium 
and also to generate a large elliptic flow vi in noncentral 
Au+Au collisions [3]. 

Since with elastic pQCD scatterings alone no ther- 
malization is expected to be achieved [2111!, the inelastic 
99 999 processes seem to play a leading role in early 
thermalization, although for the parameters chosen in [2] 
the cross section of gg — > ggg collisions is a factor of 2 
smaller than that of elastic scatterings. In order to under- 
stand this a transport cross section, 
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was introduced as a pertinent quantity measuring the con- 
tributions of different collision processes to kinetic equili- 
bration [5], since large-angle collisions should contribute 
more to momentum isotropization. The results (see Fig. 
48 in [2]) showed that even due to the almost isotropic 
distribution of the collision angles in inelastic collisions 
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its transport cross section is only the same as that of elas- 
tic scatterings. Therefore one cannot understand why the 
inclusion of the inelastic processes brings so much to ki- 
netic equilibration. At first sight it seems to be an un- 
solvable problem. On the other hand, however, there is no 
reason to believe that momentum isotropization should 
relate to the angular distribution by means of the trans- 
port cross section and not by another formula. The con- 
cept of the transport cross section may be more intuitively 
than mathematically correct. In this work we will find a 
mathematically correct way to quantify the contribution 
of different processes to thermal equilibration and com- 
pare them with each other. The core issue is the transport 
collision rate. With this quantity we "quantify" the role 
of different collision processes in kinetic equilibration. 



2 Parton cascade 



The buildup of the parton cascade is based on the stochas- 
tic interpretation of the transition rate. This guarantees 
detailed balance, which is, by contrast, difficult when using 
the geometrical concept of the cross section [5] , especially 
for multiple scatterings like ggg — > gg. The particular fea- 
ture of the numerical implementation in the parton cas- 
cade is the subdivision of space into small cell units. In 
cells the transition probabilities are evaluated for random 
sampling whether a particular scattering occurs or not. 
The smaller the cells, the more locally transitions will be 
realized. 
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The three-body gluonic interactions are described by 
the matrix element [7] 



\M gg ^ 



12. 9 2 qi 



99^999 I 



2 (ql + 
x 0(k ± A g 



il) 2 ki[(kj 
- cosh y) . 



(2) 



The suppression of the radiation of soft gluons due to the 
Landau-Pomeranchuk-Migdal (LPM) effect [gpp]. which 
is expressed via the step function in Eq. (|2|), is mod- 
eled by the consideration that the time of the emission, 
~ coshy, should be smaller than the time interval be- 
tween two scatterings or equivalently the gluon mean free 
path A g . This leads to a lower cutoff for kj_ and to a de- 
crease of the total cross section. 

In this work we simulate the time evolution of gluons 
produced in a central Au+Au collision at RHIC energy. 
The initial gluons are taken as minijets with transverse 
momentum being greater than 1.4 GeV, which are pro- 
duced via semi-hard nucleoli- nucleoli collisions. Using the 
Glauber geometry the gluon number is initially about 700 
per momentum rapidity. These gluons take about 60% of 
the total energy entered in the collision. Choosing such an 
initial condition and performing a simulation including 
Brcmsstrahlung processes we obtain dEx/dy about 640 
GeV at midrapidity at a final time of 5 fm/c, at which 
the energy density of gluons decreases to the critical value 
of 1 GeV/fm 3 . The value of dEx/dy obtained from the 
simulation is comparable with that from the experimental 
measurements at RHIC. 

We concentrate on the central region: < xt < 1.5 
fm and —0.2 < rj < 0.2, where n denotes the space-time 
rapidity. Results which will be shown below are obtained 
in this region by ensemble average. 

The importance of including inelastic pQCD gg <-> ggg 
processes to momentum isotropization is clearly demon- 
strated in Fig. [U where the time evolution of the averaged 
momentum anisotropy, < p 2 z /E 2 >, is depicted, pz and 
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Fig. 1. Momentum anisotropy. 

E are, respectively, longitudinal momentum and energy of 
a gluon. The average is computed over all gluons in the 
central region. As comparison we have also performed a 



simulation with pure elastic processes starting with the 
same initial conditions. From Fig. Q] we see that while the 
gluon system is still far from kinetic equilibrium in the 
simulation with pure elastic scatterings (dashed curve), 
the momentum anisotropy relaxes to the value at equilib- 
rium, 1/3, in the simulation including inelastic processes 
(solid curve). 

We fit the time evolution of the momentum anisotropy 
using the standard relaxation formula 



Fit) = l 



exp(- 



t-tp 
0(t ) 



(3) 



For simplicity we label now the momentum anisotropy by 
Q :=< p 2 z /E 2 >. F(t) is only equal to Q(t) at t — t . For 
fixed to the relaxation time 9 is a constant with respect 
to t. Such a fit can be done at every time point t a . The 
two thin dotted curves in Fig. Q] are fits with 9 = 0.9 fm/c 
at t = 0.3 fm/c and 9 = 2.4 fm/c at t = 1.2 fm/c. We 
find that an isotropic state is achieved at about 1.0 fm/c 
in the simulation including inelastic scattering processes. 
Moreover, we see that the relaxation time 9 is generally 
time dependent. The two values of 9 in the fits are ob- 
tained by guessing. Actually 9 can be calculated exactly, 
since in order to make a local fit one should request that 
the time derivative of F(t) and Q(t) are equal at t = to. 
This leads to 



Q(t) = F(t) 



= -(Q(*0) - Qeq.) Jhr , 
t=t V{t ) 



with Q eq = 1/3. Changing to to t gives 

Q(t) 



; eq 



-Q(t) 



l 

w 



(4) 



(5) 



Equation (O expresses the relaxation rate 1/9 of the mo- 
mentum anisotropy. In the next section we analytically 
separate the relaxation rate into different terms corre- 
sponding particle diffusion and various scattering processes, 
and we define the transport collision rate which quanti- 
fies the contribution of a certain process to momentum 
isotropization. 



3 Transport collision rate 

For evaluating the momentum anisotropy at a certain space 
point one has to go to its co-moving frame, in which we 
have 



d 3 p p 2 z 



_ ,., [2nfE 2 
Taking the derivative in time gives 



f(p,x = 0,t). (6) 



Q(t) = - 



d3p P| ^-Q(*)- 



n.l (2tt) 3 E 2 dt 



d 3 p df 
(2tt) 3 ~dt 



We replace df /dt by 

df = P 
dt E 



V/ + C 22 + C 23 + Q 



:S2 



(7) 



(8) 
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according to the Boltzmann equation. C22, C 23 and C32 
denote, respectively, the collision term of gg — ► gg, gg — > 
ggg and ggg — » gg process. It is obvious that the con- 
tribution of different processes to Q is additive. Except 
for a static system the diffusion term in Eq. ([5]) generally 
has contribution to Q(t), which we denote by Wdiff.- C 22 
has no contribution to the second integral in Eq. (|7|) due 
to particle number conservation. The same is also for the 
sum of C23 and C32 at chemical equilibrium. We rewrite 
Eq. © to 



Q(t) = W diff .(t) + W 22 {t) + W 23 (t) + W 32 (t) 



(9) 



where W 22 , W23 and W 32 correspond to the gg — ► gg, 
gg — > ggg and ggg — ► gg process, respectively. According 
to ([5|) we then obtain 



1 



6{t) 



R% } (t) + R22 it) + R 23 (*) + i?r 2 (*) , (10) 



where we define 



Qeq, ~ Q(t) 



(11) 



We see that the relaxation rate of kinetic equilibration, 
1/9, is separated into additive parts corresponding to par- 
ticle diffusion and collision processes. i?22, ^23 and R32 
stand for the transport collision rates of the respective 
interactions and quantify their contributions to kinetic 
equilibration. The extension to more than three-body pro- 
cesses is straightforward because the collision term is ad- 
ditive. We note that the definition of Rj r - in Eq. (fTT|) de- 
pends on which momentum anisotropy we are looking at. 
If one defines < \pz\/E > as the momentum anisotropy 
for instance, the form i?* r ' will change accordingly. 

Putting the explicit expression of the collision term 
via the matrix element of transition into (|8"|). we obtain 
explicit expressions for Wi, which are summarized in the 
following: 



P 2 

W 22 (t) = V, < V re l 22 > 2 -71 ( Urei-^f CT22 



W 23 (t) 



P 2 

- n < V ret 23 > 2 -71 ( V re l -^f CT 2 3 
~\ Q(t) n < v rel CT 23 >2 , 



(12) 



(13) 



W 32 (t)= -n 2 ( 732 ) , N , , , . , . 

3 ^ &E1E2E3 J 2 \ E-^ &E1E2E3 1 g 



\ n a I Viz h2 



- Q{t) n 2 



/ *» \ 

\ 8E\E 2 E 3 1 3 



(14) 



where 



x (2tt) 4 <5 (4) (pi +p 2 - p[ - p' 2 



'2' I 2 x 



(15) 



1 1 f 

~° 23 := 2^3! J dr '^ dr 2 dr ^\ M ^i'2' 3 '\ 2 x 
x (2tt) 4 <5 (4) (pi + p 2 - p'i -P2-P3) 



132 := 2! 7 ff i^^^i'2f x 

x(2^ 4 <5( 4 )(p 1+ p2+P 3 -Pi-P2) 

/ 32 -.= i / dr[dr'^\M 12Z ^ 2 ,\ 2 x 



x(2ir) 4 6^(p 1 +p 2 +p 3 -p 1 -p 2 ) 



(16) 
(17) 

(18) 



with dFi = d 3 pi/(2ir) 3 2Ei for short. 022 and 023 denote 
the standard pQCD cross section of the gg — > gg and 
99 ~ * 999 process, respectively. v re i = s/2E\E 2 is the 
relative velocity. <>a and <> 3 symbolize, respectively, 
an ensemble average over pairs and triplets of incoming 
particles. In the parton cascade simulations /(p, x, t) « 
Yli 8( 3 '(p — pi)S^(x — Xi(t)), and we can approximately 
evaluate the averages O2 and O3 in local cells which 
have small volume, but a sufficient number of (test) par- 
ticles to achieve high statistics. 

The expression of WiS in ([12]). (fl3|) and (14]) indi- 
cates the difference of the gain and loss in the momentum 
isotropization within one collision. The influence of the 
distribution of the collision angle on momentum isotropiza- 
tion is implicitly contained. However, the expression of 
the transport collision rate Rj r ' is clearly different from 
n < v re ia\ r ' >2 by the formula ([1}. (The index i denotes 
22 or 23.) Only in the special case that all particles are 
moving along the Z— axis (irrespective of ± sign) and have 
the same energy E, 

f(p,x,t) cx 6(px)S{p Y )S(pz - E) + S(p x )5(p Y )S{pz + E) , 

(19) 

the lab frame is the CM system for every colliding pair. 
In this case p 2 z /E 2 = cos 2 6* and 



R 



n < v re iaf- > 2 



(20) 



This result does not depend on the chosen direction of 
initial momentum. The only necessary conditions arc that 
all particles move along the same direction and have the 
same energy. If we define the momentum anisotropy as < 
\pz\/E >, Rf - maintains its form (f20|) . but a\ r ' is changed 
to 

. da,. 



de* 



d6* 



(i 



(21) 



Figure [2] shows the transport rates obtained from 
simulations employing the parton cascade. The solid, dashed 
and dash-dotted curves depict, respectively, the transport 
collision rates R 2 r 2 ' , R% 3 and R^ 2 ' calculated in the simula- 
tion with both elastic and inelastic collisions. We realize 
the dominance of inelastic collisions in kinetic equilibra- 
tion by computing the ratio (R 2 r 3 ' + R^/R^i ~ 5. The 
thin dotted curve in Fig. [5] presents R 22 in the simulation 
with pure elastic processes. When comparing this with the 
solid curve one cannot realize much difference. 

The contribution of particle diffusion to kinetic equi- 
libration, R%'ff , calculated in the simulation with both 
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Fig. 3. "Transport collision rate" in the concept of transport 
cross section. 



elastic and inelastic processes, is depicted in Fig. [5] by the 
symbols multiplied by —1. 



Rdiff , which is expressed by 



R 



diff- 



(*) 



1 



1 



Qeq. ~ Q(t) U 



(22) 

is not computed at a certain time as performed for the 
transport collision rates, because of the inaccurate ex- 
traction of V/. The diffusion rate is obtained by explicit 
counting of the particles which come in as well as go out of 
the central region within a time interval. Although the ex- 
traction of Rdiff has a much larger statistical fluctuation, 
the sum R%'ff, + R%2 + R^i + R32 gives a value consistent 
with the relaxation rate l/0(t) as one would calculate via 
Q{t)/(Q eq . - Q(t)) directly from Fig.[Q 

We also see that R%'ff has a negative contribution to 
momentum isotropization and is quite large. We did not 
plot —Rfaff from the simulation with pure elastic pro- 
cesses. It is slightly smaller than i? 2 2 (thin dotted curve). 
We see that there is a big difference in the diffusion rate in 
both simulations. To understand this we assume Bjorken's 
space-time picture of central ultrarelativistic heavy-ion 
collisions [TU] and use the relation derived by Baym [TTj : 



P_ —rPzdf = p z df 

E E dZ t dp z 



(23) 



Inserting Eq. (|23|) into Eq. ([22]) and calculating partial 
integrals give 



Rdiff. (t) 



-2 



(Qeq.-Q(t))t 



Q(t) 



(24) 



a comparison is necessary for understanding why the con- 
cept of transport cross section cannot explain the strong 
effect when including inelastic scattering processes. In Fig. 
[3] we depict n < v re i a tr ' >2 calculated from the cascade 
simulations. At first we look at the results in the simula- 
tion with both elastic and inelastic collisions and calculate 
the inelastic to elastic ratio (dashed versus solid curve). It 
is almost a constant around 1.7 in time, which indicates 
the dominance of the gg — ► ggg scattering processes in 
kinetic equilibration. Comparing Fig. [3] to Fig. Owe real- 
ize that the results concerning elastic scatterings have no 
strong difference. On the contrary, Rtfe is a factor of 2.5 
larger than n < v re i c^s >2- The pQCD gg — * ggg process 
is much more efficient for kinetic equilibration than one 
would expect via the transport cross section. 



4 Bremsstrahlung process and the LPM effect 

It is intuitively clear that a 2 — » 3 process will bring one 
more particle towards isotropy than a 2 — > 2 process. The 
kinematic factor should be 3/2 and appears in Eq. (flU)) 
when assuming the decompositions 



P 2 

v rei-^-cr} w Q(t) < v re i a >2 



(25) 



Analogously we compare 1^23 to W^2 . The sum of the last 
term in Eq. (fl~3|) and Eq. (fT"4|) comes from the second term 
in Eq. Q with C23 + C32 instead of df/dt and should be 
zero at chemical equilibrium: we obtain 



The formula confirms that R%'tf is always negative. Using 
the approximation < p^/E 4 >« Q 2 one can also realize 
that the larger is Q, the larger —R d r ^f is. 

We have seen that only in an extreme case the trans- 
port collision rate can be reduced to a formula directly 
proportional to the transport cross section: R tr - ~ n < 
v re i a tr - >2- It is interesting to know how the calculated 
transport collision rates differ from n < v re i a tr ' >2- Such 



Tl < V re l 023 >2 = 



1 2/ ^3 



(26) 



or equivalently i?23 
position 



3-R32 



3 \ 8E1E2E3 1 3 

Assuming further the decom 



32 



Piz 



Q(t) 



'32 



8E1E2E3 / £ 



(27) 
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we have 

3 

W 23 (t) « -(n < v re i o- 23 >2 -Q(t)n < v re ia 2 3 >i) , 

^ ^ 3 ^ 8E1E2E3 ^ ^ ^3 ^ &E1E2E3 ^ 3 

These approximate expansions together with Eq. (|26[) lead 
to VP23 « 1^32 and i?23 w 5-R32' at chemical equilibrium. 
Alone due to the kinematic reason a 2 — > 3 process is 50% 
more efficient in kinetic equilibration than a 2 — > 2 or a 
3^2 process, when 022 = 023 and 1T22 = 023 • 

The LPM effect stems from the interference of the ra- 
diated gluons (originally photons in the QED medium) 
by multiple scattering of a parton though a medium. This 
is a coherent effect which leads to suppression of radia- 
tion of gluons with certain modes (w,k). w and k denote 
energy and momentum of a gluon respectively. Heuristi- 
cally there is no suppression for gluons with a formation 
time t — wjk\, smaller than the mean free path. This 
is called the Bcthc-Hcitlcr limit, where the gluon radia- 
tions induced at a different space-time point in the course 
of the propagation of a parton can be considered as in- 
dependent events. These events within the Bethe-Heitler 
regime have been included in the parton cascade calcula- 
tions. Radiation of other gluon modes with coherent sup- 
pression completely dropped out, which is indicated by 
the 0- function in the matrix element ([2]). The inclusion 
of those radiations would speed up thermalization. How 
to implement the coherent effect into a transport model 
solving the Boltzmann equation is a challenge. 

The ©-function in the matrix element ([2|) results in 
a cut-off for fcy, the transverse momentum of a radiated 
gluon, fcx > where A g is the mean free path of a 

gluon. A higher value of the cut-off will decrease the total 
cross section of a gg — > ggg collision on the one hand and 
make the collision angles large on the other hand. The 
latter leads to a large efficiency for momentum isotropiza- 
tion. Varying the cut-off downwards to a smaller value one 
would enter into the LPM suppressed regime. Although 
parton cascade calculations set up with smaller cut-offs 
cannot completely take the LPM effect into account, one 
can roughly estimate the contribution of the coherent ef- 
fect to kinetic equilibration. Such calculations will be done 
in a subsequent paper. 



at RHIC energy the old concept of the transport cross 
sections would strongly underestimate the contribution 
of gg — > ggg to kinetic equilibration. The correct results 
showed that the inclusion of pQCD Bremsstrahlung pro- 
cesses increases the efficiency by a factor of 5 for ther- 
malization. The large efficiency stems partly from the in- 
crease of particle number in the final state of gg — > ggg 
collisions, but mainly from the almost isotropic angular 
distribution in Bremsstrahlung processes due to the ef- 
fective implementation of LPM suppression. The detailed 
understanding of the latter has to be developed in future 
investigations. 
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5 Conclusion 

Employing the parton cascade we have investigated the 
importance of including pQCD Bremsstrahlung processes 
to thermalization. The question addressed is how to under- 
stand the observed fast equilibration in theoretical terms. 
The special emphasis is put on expressing the transport 
collision rate in a correct manner. The concept of trans- 
port cross section only gives a qualitative understand- 
ing of the dominant contribution of large-angle scatter- 
ings to momentum isotropization but not a correct way to 
manifest the various contributions. In case we are study- 
ing parton thermalization in a central Au+Au collision 



